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ABSTRACT 

Wc combine surveys of the radio sky at frequencies 22 MHz to L4 GHz with data from the ARCADE- 
2 instrument at frequencies 3 to 10 GHz to characterize the frequency spectrum of diffuse synchro! rc#n 
emission in the Galaxy. The radio spectrum steepens with frequency from 22 MHz to 10 GHz. The 
projected spectral index At 23 GHz derived from the low-frequency data agrees well with independent 
iiieasureii'ieuts using only data, at frequencies 23 GHz and above. Comparing t lit- spectral index at 
23 GHz to the value from previously published analyses allows extension of t he model to higher 
Frequencies. The combined data are consistent with a power-law index 3 = — 2,64 ±9.03 at, 0.31 GHz, 
steepening by an amount A 0 = 0.07 every octave in frequency. The observed sync! not ion spectral 
steepening agrees with measurement of the cosmic ray energy spectrum ami allows extrapolation of 
synchrotron emission to frequencies 100 GHz. 

Subject headings: radio continuum: general, radiation mechanisms: non-thermal 


I. INTRODUCTION 

Synchrotron emission frpin relativistic cosmic my elec- 
trons accelerated in the Galactic magnetic field domi- 
nates the diffuse radio continuum at frequencies below 
1 GHz. It is an important foreground contaminant for 
measurements of the cosmic microwave background ra- 
diation. and also serves to probe the Galactic magnetic 
field and cosmic ray distributions. Measurements of the 
synchrotron frequency spectrum are thus of interest to 
several areas in astrophysics. 

An isotropic distribution of relativistic electrons at a 
single energy E = yme 2 propagating in a uniform mag- 
netic field B lias emissivitv 

e(v) = =-/?sinoEJ ? ’(x) , (1) 

tnc- 

where a is the pitch angle between the magnetic field 
and the line of sight,, and 

F(t) = xJ K*f 3 (xt)dxr (2) 

is defined in terms of the modi lied Bessel function of 
order 5/3 with variable .r = vjv c at id 

v c — ~ — 2 B sin o (3) 

4 tt me 

(Schwinger 1949; Westfold ct al. 1959: Cater 
1961). For a power-law distribution of electron energies 
.Y(E) \ E p propagating in a uniform magnetic field, 

the synchrotron emission is also a power law. 

T a (v)*iS (4) 


where T\ is antenna temperature, t / is the radiation fre- 
quency. and 


P-3 


2 


(5> 


(Rybicki & Llghtman 1979). 
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Measurements of the synchrotron spectral index pro- 
vide important, input for models of cosmic ray propa- 
gation. Solar modulation reduces the local cosmic my 
electron density for electron energies below a few GeV so 
that synchrotron emission provides the most direct probe 
of low-energy cosmic rays. Measurements of the’ cosmic 
ray spectrum above a few GeV in turn inform models of 
the high- frequency synchrotron spectrum. Energy losses 
from cosmic ray propagation steepen the cosmic ray spec- 
trum. increasing p toward higher energies. The observed 
steepening from p — -2.6 at 5 GeV to p ** —3.2 at 
50 GeV predicts a corresponding steepening in the syn- 
chrotron spectrum from 0 ~ “2.8 at 1 GHz to 3 ^ —3.1 
at 100 GHz (Strong, Moskalenko* Ptuskin 2007), 
Comparison of the cosmic ray spectra to the predicted 
synchrotron spectrum is complicated by confusion from 
competing radio emission sources. The diffuse radio con- 
tinuum is a super posit ion of the cosmic microwave back- 
ground, synchrotron emission, free- free emission Hum the 
warm ionized interstellar medium, and thermal emission 
from interstellar dust. Electric dipole emission from a 
•population of small spinning dust grains contributes an 
additional signal peaking at frequencies near 29 GHz 
(Draine k Lazarian 1998: AU-Hainioud er ah 2009: 
Hoang, Draine. k Lazarian 2010; Ysard k Verst mete 
2010). A number of authors have attempted to disen- 
tangle t he various emission sources to determine t he syn- 
chrotron spectral index (for a recent review see Appendix 
A of Strong, Orlando, JafFr (291 1)). Despite some dis- 
crepaut results, the general trend shows a steepening of 
the synchrotron spectrum from i ^ —2.5 at 22 MHz to 
3 — - 3.0 above 23 GHz, in rough agreement with the 
observed cosmic my spectra. 

Several factors contribute to the observed scatter in 
estimates of the synchrotron spectral index. Most es- 
timates, particularly those below 23 GHz, assume a 
power-law spectrum for synchrot ron and do not explic- 
itly model spectral steepening. Comparisons between 
closely- separated frequencies more arcurafctdy reflect ilic 
local synchrotron spectrum, but have larger uncertain- 
ties from competing emission sources or measurement 
offsets. Analyses with broader frequency coverage refiner 
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foreground and offset imccrtainties but average over any 
spectral sl eepeuiiig. 

Two additional effects are important fen- analyses in- 
el tiding data from the Wilkinson Microwave Anisotropy 
Probe (WMAP) at frequencies 23 to 94 GHz. A grow- 
ing body of evidence suggests t hat a substantial fraction 
of the diffuse continuum near 23 GHz consists of electric 
dipole radiation from a population of small, rapidly spin- 
ning dust grains (Kogut et al. 199(5; dc Oliveira-Cost a 
et al. 1997. 2004; Mi ville-Deschenes et al. 2008: Dobler 
k FI uk twiner 2008; Ysanb Mi ville- Deschenes . k Ver- 
st raete 201(1: Kogut et ah 2011; Gold et ah 2011: 
Planck collaboration 2011). Analyses that, ignore this 
component tend t o over- predict the synchrotron ampli- 
tude near 23 GHz, biasing the derived spectral index to 
flatter values when comparing to lower frequencies and 
steeper values when comparing to higher frequencies. 

A second systematic error can result from improper 
treatment of offsets in the data. Measurements from 
radio surveys at frequencies below 20 GHz include tlie 
absolute intensity (zero level) of the sky. The WMAP 
differential radiometers arc insensitive to any constant 
(monopule) intensity on the sky: the zero level of the 
WMAP sky maps is set so that the map intensity in the 
Galactic polar caps matches a cosecant fit to the mid- 
latitude sky (Bennett et ah 2003: Hinshaw et al. 2009)* 
Analyses that directly compare low- frequency radio sur- 
veys to the WMAP data without subtracting a monopule 
component from the radio data will miss the correspond- 
ing emission hi the WMAP bands, biasing the derived 
spectral index to steeper values. 

The zero level bias can be significant* The H aslant 
et ah (1981) survey at 408 MHz is commonly used to 
model synchrotron emission. The North Galactic pole 
has measured temperature 19 ±3 K at 408 MHz, while 
a esc \ h\ fit. to the same 408 MHz map predicts a polar 
contribution of only 5.1 ± 0.6 K. Similar results apply t o 
t he Southern hemisphere, where the measured polar cap 
temperature of 21 ± 3 K significantly exceeds the value 
4.0 ±0.5 K obtained from a csc|5| fit. Only 2.7 K of the 
difference can be attributed to emission from the eosmie 
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FlU. 1 Toy model thawing effect of improper zero I eve] *ub- 
irat’tioii tm the derived .spectral index between 408 MH* and 23 
GHz. 


TABLE 1 

Sky Suhvkyf Used fob Synchrotron Analysis 


frequency 

(Gil*) 

Calibration 
U ueerlainty 

Offset 

Uncertainty (K) 

Relative 

Uncertainty" 

0.022 

n.u'j 

3000' 

0.15 

0.045 

O.10 

250 

0.11 

0.408 

0.10 

3.0 

0.17 

1.420 

0.05 

0.5 

0.63 

3/20 

0.001 

0.011 

0.10 

3,41 

0.001 

0.006 

0.07 

7.98 

u.om 

0,036 

0.89 

8.33 

0,001 

0.042 

2.64 

9.72 

0.001 

0,003 

0.34 

10,49 

0.001 

0.002 

0.27 


^Quadrature sum of calibration and offset uncertainties, divided 
by the mid-latitude sky temperature. 

microwave background, leaving a large residual. 

Figure 1 illustrates the bins induced by including this 
residual at 408 MHz but excluding it From the WMAP 
data. Wo take the 408 MHz map, remove the 2.7 K 
CMB monopule, and scale the remaining radio emission 
to 23 GHz using a power- law index 3 = -2.7, We then 
remove a mo no pole from the scaled map so that the map 
temperature in the south polar cap matches the cse|5| 
fir. We compute the bias in apparent spectral index In- 
comparing the 408 MHz map to the scaled 23 GHz map 
before and after removing the seated monopole. Figure 
1 shows the bias in spectral index a] 23 GHz. binned 
by Galactic latitude. Dis-simikir treatment of the map 
zero level creates a spatially varying bias Ad m 0,15, 
comparable to the total spectral steepening predicted by 
the measured cosmic ray spectra. The bias is largest 
in regions where the sky brightness is faintest, at high 
latitudes or away from the Galactic center. 

Measurement uncert.tunties in the absolute level of the 
sky brightness can also introduce bias in estimates of the 
synchrotron spectral index. Many of the low-frequency 
surveys have uncertainty in the measured zero level ap- 
proaching 30% of the polar cap brightness. As with 
the toy model above, sucli measurement errors intro- 
duce spatially dependent, biases that are largest where 
the sky brightness is faintest. Minimizing uncertainty in 
the derived synchrotron spectrum requires a combination 
of measurements with good sky coverage and good con- 
trol of offset uncertainty at frequencies where competing 
emission sources are faint. No such ideal data set yet 
exists. In this paper, we model the synchrotron spectral 
index and curvature using low-frequency radio surveys 
with high sky coverage but large offset uncertainty, com- 
bined with higher- frequency measurements with limited 
sky coverage but excellent offset to certainty. 

2. SKY MARS 

We model synchrotron emission using radio surveys at 
22 MHz (Roger et al, 1999), 45 MHz (M aeda et al. 1999: 
Alvarez et al. 1997), 498 MHz (Haslam et al. 1981). and 
1420 MHz (Reich, IcstorL k Reich 2001; Reich k Reich 
1986). These surveys have full or nearly- full sky cover- 
age at frequencies where Galactic radio emission is sig- 
nificant, with gain and zero-level systematica controlled 
at the 19 20% level. We supplement t he radio surveys 
with sky maps from the Absolute Radiometer for Cos- 
mology. Astrophysics, and Diffuse Emission (ARCADE 
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FIG. 2, Sky temperature* and best-fit model (dashed line) 
for a 4 11 diameter patc h on the Galactic plane centered on the 
the brightest pixels in the ARCADe 2 sky coverage. For clarity, 
the data are plotted relative to the hest-fii power-law model (d = 

2.56) with zero curvature (solid line). All tits include the non- 
trivial covariance between individual data points. The data are 
well descr ibed by a model with spectral index 3 - -2.60 ± 0.04 
and curvature C - -0.081 ±0.028. 

2) instrument 2 at 3. 8, anti 10 G Hz (Kogut H al. 2011). 
Tlic ARCADE 2 data observe both the Galactic plane 
and mid-latitude regions (|6j < 10°) with sufficient con- 
trol of zero-level uncertainty to constrain the synchrotron 
curvature relative to the lower- frequency radio surveys. 

Table 1 sun tin arises the input sky maps, The increase 
in the offset uncertainty at low frequency is compensated 
by a corresponding increase in sky brightness* The fi- 
nal column shows the relative measurement uncertainty 
for a mid-latitude region, defined as the ratio of the 
combined offset and calibration uncertainty to the mea- 
sured brightness at (/.ft) — {17° f — 35°) after removing 
the CMB monopole. The selected maps provide roughly 
uniform relative sensitivity to synchrotron emission over 
2,5 decades of frequency. 

VVe convert all maps to units of antenna temperat ure 
mid subtract the CMB monopole at ( thermodynamic ) 
temperature 2.725 K from the measured sky tempera- 
tures. We then convolve each map to the 11? 6 angular 
resolution of the ARCADE 2 instrument. At frequen- 
cies of 10 GHz and below, both thermal dust emission 
and spinning dust emission are negligible. Free-free emis- 
sion, however, can still be appreciable. We correct the 
convolved maps by scaling the WMAP 7-year maximum 
entropy free-free model (Gold rt aL 2011) to each Fre- 
quency using spectral index —2.15, convolving t he scaled 
model to ll?fi angular resolution, and subtracting the 
convolved model from each sky survey. The resulting 
maps arc- dominated by synchrotron emission. 


3. ANALYSIS 

The input sky maps define a data set T(ili>) sampled 
at discrete pixel directions/i and 10 discrete frequencies v 
ranging from 22 MHz to 10 GHz, We model synchrotron 

2 The ARCADE data are available at the Legacy 
Archive For Microwave Background Data Analysis. 
http://lambda.gflf c .naBa.gov 


emission as a modified power law 

/ v 

( 6 ) 

with sj H'vt ral n u*va l 1 1 re C dt *fi i ied wit h res port tc * n Ter- 
ence frequency I'n - 310 MHz. The adopted value for 
vu minimizes covariance between the fitted amplitude .1 
and spectral index ft* simplifying extrapolation to other 
frequencies. 

For each pixel h we define a 10 by 10 data covariance 
matrix M with diagonal elements determined by rhe in- 
strument noise, calibration mid offset uncertainty. The 
input maps are not all linearly independent. Measure- 
ments at 22 MHz used the 408 MHz map to determine 
the declination dependence of the gain. YVe estimate the 
result ing correlation of spatial structure in the two maps 
at 50%, The ARCADE 2 maps have independent instru- 
ment noise but share a fraction of the offset uncertainty 
related to absolute themometry uncertainty and ground 
glint (Singal et ai. 2011). All 10 maps share a common 
model for free- free emission. We conservatively estimate 
the uncertainty in the free-free correction at 30% of the 
free- free amplitude; however, the results do not change 
significant ly as the model free-free amplitude is varied 
by as much as 50% * Off-diagonal elements in M include 
these effects. 

We restrict the analysis to the 8% of the sky observed 
at all 10 frequencies. For each pixel, a least-squares 
minimization determines the best-fit parameters A , J, 
and C, Figure 2 shows the measured temperature and 
best-fit model for the brightest Galactic plane region 
[Lb) (— 52°,0*) within the common sky coverage. The 

t \ ata slit i w e v i d enc c U >r spt ( t ra 1 c i i r vat ure , with best- fi I 
values ft = -2.60 ± 0,04 and C - -0,081 ± 0.028 evalu- 
ated at uq 310 MHz. The spectral curvature C for this 
region is significant at approximately 3 standard devia- 
tions compared to the baseline model with C — l), 

We may extrapolate the spectral models to compare 
the results at frequencies II) GHz and below to indepen- 
dent determi nations of the spectral index using WMAP 
data at frequencies 23 GHz and above. We compute 
the antenna temperature of the modeled spectra to de- 
rive the effective power- law index for frequencies near 23 



FIG* 3.— Spectral index evaluated at 310 MHz from the 11* 
Frequency fit . 
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I r I<:. 4. Spectral curvature from U -frequency fit. The mean 
curvature C — —0.052 corresponds to a steepening of the locaf 
spectrwi index bv an amount Af3 = 0.07 every octave in frequency. 

GHz* Note that this is nut equivalent to evaluating Eq. 
G at w 23 GHz. whirl i would yield the sealing from 310 

MHz to 23 GHz but not the power-law index at 23 GHz. 
The mean for all 258 pixels in the common sky coverage 
is i%% = -3.02 with standard deviation 0.22* 

The extrapolated value com pares well with indepen- 
dent determinations of the spectral index above 23 GHz. 
Kogut et ah (2004) analyze WMAP polarization data 
to derive synchrotron spectral index 3 = —3.2 ± 0*1 av- 
eraged over the hill sky, Dunkley et al. (2009) use a 
Bayesian analysis of polarization data and find the mean 
spect ral index 3 = -3.03 ±0.01 with pixel- to- pixel stan- 
dard deviation 0.25 the high-latitude sky. Gold et ah 
(2011) use template fitting techniques to derive spectral 
index 3 = —3.13 between 23 and 33 GHz. 

Much of the scatter in t he extrapolated spectral indices 
results from pixels at high latitude where the emission is 
faintest. The extrapolated index in these pixels can reach 
i inphysical values. We reduce the scat t er by constrain- 
ing the spectral index at 23 GHz. For each pixel, we use 
the radio data (Table l ) to fit the synchrotron amplitude 
A(tt) over a 2-dimensional grid in t lie spectral paramet ers 
J and C. At each grid point, we compute the \~ value 
/? 7 A / “ 1 7? where M~ l is the inverse covariance matrix 
and /? is the difference vector between the measured and 
i node Jed temperatures. We then use die spectral param- 
eters 3 and C to model the power- law index at 23 GH/ 
and compare the resulting value to a prior. We use the 
difference between the extrapolated spectral index and 
the prior to augment tile \ 2 at each grid point. 



where 3%:\ is the model spectral index evaluated at 23 
GHz, and 3 t , ±a v = - 3.1 ± 0.1 is the prior at 23 GHz. 
The minimum \ 2 over the entire grid then defines the 
best-fit model at that pixel* This allows inclusion of the 
spectral information derived from frequencies above 23 
GHz without confusion from either additional emission 
components {spinning dust ) above 23 GHz or the missing 
zero level in the WMAP dat a. 

Figures 3 and 4 show t he distribution of the resulting 


best- fit values for the api'ctm! index and curvature* In- 
cluding the constraint at 23 GHz, the best-fit spectral 
index has mean value 2.64 and standard deviation 0*03 
f 1 1 rc iference free [ ue 1 1 cy t/q 310 MHz. The 1 ) est-fi t cur- 

vature has mean value -0.052 with standard deviation 
0.005* The corresponding spt'ctral index at 23 GHz is 
(fas) = -3.09 with standard deviation 0.05. Inclusion of 
the prior at 23 GHz does not induce a significant shift in 
the mean for the extrapolated spectral index, but does 
significantly reduce the pixel- to- pixel scatter, 

4. DISCUSSION 

Radio data show statistically significant steepening of 
the synchrotron spectrum from 22 MHz to It) GHz. The 
nearly uniform relative uncertainty of the selected data 
minimizes dependence of the fitted parameters cm off- 
set or calibration errors at any one frequency. We test 
whether the best- fit parameters are particularly sensi- 
tive to any one input map by repeating the analysis af- 
ter dropping one or two maps from the fit. We select 
either one radio survey (22 MHz. 45 MHz* 408 MHz, 
or 1420 MHz) or a pair of ARCADE frequency chan- 
nels (3 GHz* 8 GHz* or 10 GHz) and repeat the fit af- 
ter deleting the corresponding elements from the data 
vector T and covariance mat rix AL The resulting shift 
in either the spectral index or curvature parameters is 
smaller than the pixel-to-pixel standard deviation using 
all 111 frequency channels. Systematic enors in the offset 
or temperature calibration do not appear to dominate 
the mill ti- frequency analysis. 

The steepening of the synchrotron spectrum is con- 
sistent with models of cosmic ray propagation in the 
Galactic magnetic field, Jaffe et ah (2011) combine ra- 
dio observat ions wit h the gauprop 3 cosmic ray propaga- 
tion code to model synchrotron emission on the Galactic 
plane. They find a power-laiv index -2.8 < 3 < -2.7 J 
from 408 MHz to 2.3 GHz and -2.98 < 3< -2.91 from 
2.3 GHz to 23 GHz, The corresponding values for the 
spectral steepening model are 3 = 2.7fi from 408 MHz 

to 2.3 GHz and 3 -2,97 from 2.3 GHz to 23 GHz. 

We may use the best-fit values in each pixel to predict 
l he synchrotron spectrum at higher frequencies where t he 
emission is fainter and competing sources stronger. Pre- 
vious attempts to disentangle competing emission from 
free- free, synchrotron, thermal dust, and spinning dust 
emission using WMAP data have suffered from degen- 
eracy between the synchrotron and spinning dust emis- 
sion* both of which are falling at frequencies above 33 
GHz (see. c.g.. the discussion in Gold et ah (2011))* 
Extending the synchrotron curvature observed at lower 
frequencies into the millimeter baud reduces confusion 
bet w *et i the sj mini ug d us t and sy n cl in >i rc n i spt -c -t ra ant 1 
may facilitate diarctemat-ion of both the spatial dist ri- 
bution and frequency sptvt rum of spinning dust emission 
in the interstellar medium* 

Table 2 shows the local j lower- law index T x v 1 at se- 
lected frequency bands. The modeled spectrum steepens 
by A3 = 0*07 every octave in frequency. from 3 = 2.67 

at 408 MHz to 3 — -3.24 at 94 GHz* Note* however* 
that the spectral steepening observed at low frequencies 
can not continue indefinitely, Ftrmi measurements of 
the cosmic ray energy spectrum are consistent with a 

:i http ;//galprop , Stanford . edu 
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TARLE 2 

Local Power-Law Spectral Inuex 


Frequency 

(mz) 

Power-Law 

Index 

0.408 

-2.67 

23 

-3.09 

:i3 

3.13 

4i 

-3.15 

61 

-3.39 

94 

-3.24 


single power law from energy 7 GeV to 1 TcV (Abdo <?T 
al* 2009; Ackerniami et al. 2010, 2012)* If anything, the? 
Fermi data suggest a .modest fiat truing of the cosmic ray 
energy spectrum at higher energies. which would induce 1 
a positive curvature to the synchrotron spectrum. Direct 
confirmation of the synchrotron spectrum above 23 GHz 
remains a challenge* 

5, CONCLUSIONS 

Radio data are consistent with a synchrotron spectrum 
that steepens with frequency from 22 MHz to 10 GHz. 
Direct comparison of low-frequency radio surveys with 


the WMAP data at 23 to 94 GHz is complicated both 
by the presence of additional emission components at 
higher frequencies and by the subtraction of a substan- 
tia) monopole component of sky emission by the differen- 
tial WMAP instrument. The synchrotron spectral index 
at 23 GHz* derived using only lower-frequency radio sur- 
veys. is consistent with the value derived independently 
using only dat a at higher frequencies. We extend t he ra- 
dio data by comparing the extrapolated index at 23 GHz 
to a prior based on higher- frequency dat a* The combined 
data have mean spectral index 3 -2*64±0.03 and cur- 

vature C = -0.052 ± 0.005 at. reference frequency 0,31 
GHz. The measured spectrum steepens by an amount 
Aii 0.07 every octave in frequency. The observed syn- 
chrotron spectra! steepening agrees with measurement of 
t he cosmic ray energy spect rum and allows extrapolation 
of synchrotron emission to frequencies 100 GHz. 

This research hs based upon work supported by the 
National Aeronautics and Space Admiiiistmt ion through 
the Science Mission Directorate under the Astronomy 
and Physics Research and Analysis suborbit al program. 
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